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(57) Abstract 

PURPOSE: To efficiently calculate an optimal design 
value by properly setting variable difference amount in 
the design of lens, etc., by the method of descent 

CONSTITUTION: A CPU 1 calculates the eigen value 
distribution information of the product matrix (AvTAv) 
between the Jacobfan matrix Av having the variable 
amount of the characteristic (evaluation function Fv(Xv)) 
for variable difference amount Sxj of the system being 
the design object as an element and transpose of 
matrix AvT and sets the variable difference amount SXj by 
using the eigen value distribution information. Instead 
of the eigen value of the product matrix (AvTAv), the 
square value of the singular value of the Jacobtan 
matrix A may be used. 
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8 (1/n) = 1x10- ,5 (1/nO =1X10"* 
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8 (1/n) = 1x10-, 5 (l/i2) =1X10° 

si = 5.5 x 10- ,sz = 6.8 x 10° (S3 (d)) 
8 (1/n) =1X10"- ,5 (I/12) = 1x10- 

si - 8.6 x lO' 18 ,S2 = 2.2 x 10- (04 (a)) 
8 (1/n) = lxlO- 4 .5 (1/nO = 1 x 10" 4 

si = 8.6 x 10- t sz = 2.2 x 10-" (0 4 (b)) 
5 (1/n) =1x10-, 5 (I/12) = 1x10- 

si = 8.2 x 10- ,S2 - 2.3 x 10" a (0 4 (c)) 
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si = 5.5 x 10'" ,sz = 6.8 x 10° (04 (d)) 
5 (1/n) = 1x10-.$ (l/r2) = 1x10- 

si = 8.6 X 10" 18 ,S2 = 2.2 x 10° (0 5 (a)) 
8 (1/n) = 1x10-, 5 (1/X2) =1X10-* 

si = 8.6 x 10-" ,S2 = 2.2 x 10° <M 5 (b)) 
8 (1/n) = lxlO-,<5 (1/nO = 1x10- 

si = 2.2 x 10° .62 = 8.7 x 10" B (0 5 (c)) 
8 (1/n) = lxlO-,<5 (l/r2) =1X10° 
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8 (1/n) =1x10° ,8 (1/nO = 1x10- 
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si = 5.9 x 10- ,S2 — 5.2 x 10- (0 6 (c) 
8 (1/n) = lxlO°,£ (1/rz) =1x10° 

si = 1.5 x 10° .52 = 5.2 x 10- (B 6 (d)) 
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do= 1 • 0 (2 3) 

[0 0 3 9] 
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d B XlO" B ~d,xO. 5X10* (24) 
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5 (1/n) = 1x10-. <5 (1/n) = 1x10- (B8 (c)) 
8 (1/n) = lxl0-,<5 (1/rz) =1x10° (03 (d)) 
5 (1/n) = lxi0-,<5 (1/ia) = 1x10- (04 (a)) 
8 (1/n) = lxl0-,<5 (1/rz) =1x10- (04 (b)) 
8 (1/n) = 1x10-, 5 (1/n) =1x10- (05 (a)) 
8 (1/n) =lxl0- 3 .<5 (1/rz) = 1 X 10- (B5 (b)) 
5 (1/n) =1x10°. 5 (l/n) =1X10- (06 (a)) 
5 Q/ n ) =1x10° .8 (1/rz) = 1X10- (06 (b)) 

(26) 

SfcB Wffl s 1 , s 2 tlTO * Q jMUT 
( (2 7) it) T*£><5o 
[0 0 4 2] 
[82 3] 

8 (1/n) =1x10-, 8 (1/rz) = 1 x 10- (04(c)) 
5 Q/ ri ) = ixl0-,<5 (1/rz) =1X10° (04 (d)) 
$ <i/ n ) = 1X10'* ,6 (1/rz) =1X10- (05 (c)) 
8 (1/n) = lxl0-,£ (1/rz) =1x10° (05 (d)) 
5 (l/„) = l x 10° ,5 (1/rz) = 1 x 10- (06 (c)) 
8 (1/n) = lxl0°.<5 (1/rz) =1x10° (0 6(d)) 
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(27) scicsufct^n^oK»«»**iafi'3r«c 

<h ICJ: U tt«SttMIE3Rd O«*»*r*fci0OKaft 

[0 0 4 3] «ky W5Slcfl9Ea8M»«(Oii^ffi*ia 
Kr*icii, mettttJEadOttmttdotf. *frM 

(A V U V ) OB 3 Stt# : «i©«^tOt)OJ:U'hS<«/h 
cD ! t)<^cfcy^^<^:oTU^t)(7)^iS:^:-rn(Sci:iN. * 

[0 0 4 4] 
[82 4) 

$ (i/ n ) = ixl0-, <5 (1/rz) =1x10° (03 (d)) 

$ d/n) =1x10-, 8 (1/rz) =1x10° (04 (d)) 

5 (1/n) =1X10-, <5 (1/rz) = 1X10- (05 (a)) 

8 (1/n) = lxl0-,£ (1/rz) = 1x10- (05 (b)) 

8 (1/n) = lxl0-,<5 (1/rz) =1x10° (06 (c)) 
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[0 0 6 0] 
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6 ( B ) =ixio- a -ixi 0° (31) 

[0 0 6 1 ] 
[tt2 8] 

7" x 4 11 = 8.3 x 10* A & 

(Genetic Algorithm) ^S£*«fcU»6nxi^, it£ 
tt7;U=f>JXAI41 9 6 O^ttlCJ. KHollandlCcfcU^ 

[0 0 6 2] AteMlctem 2 icsrj:5*¥JII"C«7 
$ti^o tf^f^^S 1 TttWBt^E^WkOTSie 

(loiooiio) 3ttt # *flE5. Jfcic, gfc^ite 1 ? 1 
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ISJBttWIEaotafflltcloOTSIfTM (A V TA V ) OBWfB 
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[0 0 6 7] J:UIKBltiif5tt*IMIiiO€k5i-r««* 

tt— Kl^lMIS^»6ti*J:5Cr*. 
[0 0 6 8] **JE«Ttt. 01 0IC*-r«fc5ft#ffi© 
U IMEB*B#d«»«B*W:U'b«B 

[0 0 6 9] AftttttttttflTKH-rfia:. **TBMt 

r £®BPBg t**tl**l*illC 

[0 0 7 0] 
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6 (1/r) =1X10' 6 

a (g) =1X10* 3 (33) 



If 1 0 1 1 0 0 0 1, 01101010, 111011 
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(KM) *W£©lWBtt*fiH*TfM»U VXDXtv 
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3-3-3-3-3-3-3-3 (32) 

£fc<o HMEB^rtWM&tB^*^ SI 1 

(a) lc5Vr«fc5lC£:o*:e CO*lttt*JftiLT» 
E!11 (a) OJ:5«:B*tt»*lcH1 O^cfcp^K*- 

gr, ®F e 1PSglc•^t^T^■n J enlUT(Oct^lC^^^lo 
[0 0 7 1 ] 

[ft 3 0] 

6 (1/r) =1 XI 0" B ~1X1 0" 2 

5 («) =lX10" a ~lXl0" 1 (34) 

C©«»8#lK£tl (3 4) lcS*j£-r*B*<t#fii«: 
WKStHi i (b) ic^-TckPiefco*:. Hi i 
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6 (1/n) = 1x10-, 5 (l/rz) = 1x10- 

Pl 2 = 8.6 x 10~» ,p* 2 = 2.2 x ID" 8 3 (a)) 
6 (1/n) = lxlQ-,5 (1/nO = 1x10- 

p, 2 = 8.0x 10- ,0 ,pa 2 = 2.3x10- (03 (b)) 
6 (1/n) = lxl0-,S Q/r2) - 1 x 10- 

P! 2 = 9.0x 10" 10 ,P2 2 =2.1 x 10"" (03 (c)) 
6 (1/n) = 1x10-,$ (J/rz) =1x10° 

Pl 2 = 5.5 x 10"* fPa 2 = 6.8 x 10° (S 3 (d)) 
6 (1/n) = 1x10-, 5 C1/T2) = 1x10- 

p, 2 «8.6 x 10- 18 tPa 2 = 2.2x 10- (04 (a)) 
5 (1/n) =1x10-, 6 (1/nD = 1x10- 

p, 3 = 8.6 x 10- B .p, 2 = 2.2 x 10- (S14 (b)) 

5 (1/n) = 1x10- ,5 G/i2) = 1x10- 

Pl 2 = 8.2x 10- ^=2.3 x 10- (04 (c)) 

6 (1/n) = lxlO-,<5 G/nO =1x10° 
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(54) [Title of the Invention] METHOD AND DEVICE FOR 
SUPPORTING DESIGN OF LENS OR THE LIKE 

(57) [Abstract] 

[Object] To efficiently calculate an optimal design value 
by adequately setting a variable difference amount in the 
design of a lens or the like by an attenuation least square 
method. 

[Solving Means] A CPU 1 calculates eigen value 
distribution information on the product matrix (Av T Av) of the 
Jacobian matrix Av with the variable amount of the 
characteristic (an evaluation function Fv(Xv) ) of the system 
for the variable difference amount 6xj of the system being a 
design object as elements thereof and the transposed matrix 
Av T and sets the variable difference amount 6xj by using the 
eigen value distribution information. The square value of 
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the singular value of the Jacobian matrix A may be used in 
place of the eigen value of the product matrix (Av T Av) . 
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[Claims] 

[Claim 1] A design supporting device to support the design 
of a lens or the like by a descent method for locally- 
changing and optimizing a plurality of components so that a 
system of a design object comprising the plurality of 
components has a desired characteristic comprising: 

an eigen value distribution information operating means 
to operate eigen value distribution information on the 
product matrix (Av T Av) constituted of the product of the 
Jacobian matrix (Av) with the variable amount of the 
characteristic of the system to the difference amount of the 
variables corresponding to the plurality of components as 
the elements thereof and the transposed matrix (Av T ) of the 
Jacobian matrix; and 

a difference amount setting means to set the difference 
amount of the variables by using the eigen value 
distribution information. 

[Claim 2] The design supporting device according to Claim 1, 
wherein the difference amount of the variables is set so 
that the initial value of the non-linear correction term 
corresponding to the variable for an optimizing objective 
function is present between the minimum value and the 
maximum value of the eigen value distribution of the product 
matrix. 

[Claim 3] The design supporting device according to Claim 1, 
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wherein the difference amount of the variables is set so 
that the median of the eigen value distribution of the 
product matrix is substantially matched with the median of 
the initial value of the non-linear correction term 
corresponding to the variables for an optimizing objective 
function. 

[Claim 4] The design supporting device according to Claim 1, 
wherein the difference amount of the variables is set so 
that the mean of the eigen value distribution of the product 
matrix is substantially matched with the median of the 
initial value of the non-linear correction term 
corresponding to the variables for an optimizing objective 
function. 

[Claim 5] The design supporting device according to Claim 1> 
wherein the difference amount of the variables is set so 
that the median of the eigen value distribution of the 
product matrix is substantially matched with the mean of the 
initial value of the non- linear correction term 
corresponding to the variables for an optimizing objective 
function. 

[Claim 6] The design supporting device according to Claim 1, 
wherein the difference amount of the variables is set so 
that the mean of the eigen value distribution of the product 
matrix is substantially matched with the mean of the initial 
value of the non- linear correction term corresponding to the 
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variables for an optimizing objective function. 
[Claim 7] The design supporting device according to Claim 1, 
wherein the difference amount of the variables is set by- 
using eigen value distribution information on the product 
matrix and predetermined weighting information set 
corresponding to the eigen value distribution. 
[Claim 8] A design supporting device to support the design 
of a lens or the like by a descent method for locally 
changing and optimizing a plurality of components so that a 
system of a design object comprising the plurality of 
components has a desired characteristic comprising: 

a square value distribution information operating means 
to operate square value distribution information on the 
eigen value of the Jacobian matrix (Av) with the variable 
amount of the characteristic of the system to the difference 
amount of the variables corresponding to the plurality of 
components as the elements thereof; and 

a difference amount setting means to set the difference 
amount of the variables by using the square value 
distribution information. 

[Claim 9] The design supporting device according to Claim 8, 
wherein the difference amount of the variables is set so 
that the initial value of the non-linear correction term 
corresponding to the variables to an optimizing objective 
function is present between the minimum value and the 
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maximum value of the square value distribution of the eigen 
values of the Jacobian matrix. 

[Claim 10] The design supporting device according to Claim 
8, wherein the difference amount of the variables is set so 
that the median of the square value distribution of the 
eigen value of the Jacobian matrix is substantially matched 
with the median of the initial value of the non- linear 
correction term corresponding to the variables to an 
optimizing objective function. 

[Claim 11] The design supporting device according to Claim 
8 # wherein the difference amount of the variables is set so 
that the mean of the square value distribution of the eigen 
value of the Jacobian matrix is substantially matched with 
the median of the initial value of the non-linear correction 
term corresponding to the variables to an optimizing 
objective function. 

[Claim 12] The design supporting device according to Claim 
8, wherein the difference amount of the variables is set so 
that the median of the square value distribution of the 
eigen value of the Jacobian matrix is substantially matched 
with the mean of the initial value of the non- linear 
correction term corresponding to the variables to an 
optimizing objective function. 

[Claim 13] The design supporting device according to Claim 
8, wherein the difference amount of the variables is set so 
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that the mean of the square value distribution of the eigen 
value of the Jacobian matrix is substantially matched with 
the mean of the initial value of the non-linear correction 
term corresponding to the variables to an optimizing 
objective function. 

[Claim 14] The design supporting device according to Claim 
8, wherein the difference amount of the variables is set by 
using square value distribution information on the eigen 
value of the Jacobian matrix and predetermined weighting 
information set corresponding to the square value 
distribution . 

[Claim 15] A design supporting method to support the design 
of a lens or the like by a descent method for locally 
changing and optimizing a plurality of components so that a 
system of a design object comprising the plurality of 
components has a desired characteristic comprising: 

a step of operating eigen value distribution 
information on the product matrix ( A v T A v ) constituted of the 
product of the Jacobian matrix (Av) with the variable amount 
of the characteristic of the system to the difference amount 
of the variables corresponding to the plurality of 
components as the elements thereof and the transposed matrix 
(Av T ) of the Jacobian matrix; and 

a step of setting the difference amount of the 
variables by using the eigen value distribution information. 
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[Claim 16] The design supporting method according to Claim 
15, wherein the difference amount of the variables is set so 
that the initial value of the of the non-linear correction 
term corresponding to the variables to an optimizing 
objective function is present between the minimum value and 
the maximum value of the eigen value distribution of the 
product matrix . 

[Claim 17] The design supporting method according to Claim 
15, wherein the difference amount of the variables is set so 
that the median of the eigen value distribution of the 
product matrix is substantially matched with the median of 
the initial value of the non-linear correction term 
corresponding to the variables to an optimizing objective 
function . 

[Claim 18] The design supporting method according to Claim 
15, wherein the difference amount of the variables is set so 
that the mean of the eigen value distribution of the product 
matrix is substantially matched with the median of the 
initial value of the non-linear correction term 
corresponding to the variables to an optimizing objective 
function. 

[Claim 19] The design supporting method according to Claim 
15, wherein the difference amount of the variables is set so 
that the median of the eigen value distribution of the 
product matrix is substantially matched with the mean of the 
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initial value of the non-linear correction term 
corresponding to the variables to an optimizing objective 
function . 

[Claim 20] The design supporting method according to Claim 
15 , wherein the difference amount of the variables is set so 
that the mean of the eigen value distribution of the product 
matrix is substantially matched with the mean of the initial 
value of the non-linear correction term corresponding to the 
variables to an optimizing objective function. 
[Claim 21] The design supporting method according to Claim 
15, wherein the difference amount of the variables is set by 
using eigen value distribution information on the product 
matrix and predetermined weighting information set 
corresponding to the eigen value distribution. 
[Claim 22] A design supporting method to support the design 
of a lens or the like by a descent method for locally 
changing and optimizing a plurality of components so that a 
system of a design object comprising the plurality of 
components has a desired characteristic comprising: 

a step of operating square value distribution 
information on the eigen value of the Jacobian matrix (Av) 
with the variable amount of the characteristic of the system 
to the difference amount of the variables corresponding to 
the plurality of components as the elements thereof; and 

a step of setting the difference amount of the 
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variables by using the square value distribution information. 
[Claim 23] The design supporting method according to Claim 
22, wherein the difference amount of the variables is set so 
that the initial value of the non- linear correction term 
corresponding to the variables to an optimizing objective 
function is present between the minimum value and the 
maximum value of the square value distribution of the eigen 
values of the Jacobian matrix. 

[Claim 24] The design supporting method according to Claim 
22, wherein the difference amount of the variables is set so 
that the median of the square value distribution of the 
eigen value of the Jacobian matrix is substantially matched 
with the median of the initial value of the non- linear 
correction term corresponding to the variables to an 
optimizing objective function. 

[Claim 25] The design supporting method according to Claim 
22 , wherein the difference amount of the variables is set so 
that the mean of the square value distribution of the eigen 
value of the Jacobian matrix is substantially matched with 
the median of the initial value of the non-linear correction 
term corresponding to the variables to an optimizing 
objective function. 

[Claim 26] The design supporting method according to Claim 
22, wherein the difference amount of the variables is set so 
that the median of the square value distribution of the 
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eigen value of the Jacobian matrix is substantially matched 
with the mean of the initial value of the non- linear 
correction term corresponding to the variables to an 
optimizing objective function. 

[Claim 27] The design supporting method according to Claim 
22, wherein the difference amount of the variables is set so 
that the mean of the square value distribution of the eigen 
value of the Jacobian matrix is substantially matched with 
the mean of the initial value of the non-linear correction 
term corresponding to the variables to an optimizing 
objective function . 

[Claim 28] The design supporting method according to Claim 
22, wherein the difference amount of the variables is set by 
using square value distribution information on the eigen 
value of the Jacobian matrix and predetermined weighting 
information set corresponding to the square value 
distribution. 

Detailed Description of the Invention 
[0001] 

[Technical Field of the Invention] The present invention 
relates to a design supporting device and a design 
supporting method using an attenuation least square method 
out of optimization methods to locally changing a plurality 
of components so that an objective system comprising the 
plurality of components have desired characteristics, and 



- 12 - 



for example, it relates to a design supporting device and a 

design supporting method to support a lens design, an 

optimum pattern design of an LSI, a building design and the 

like. 
[0002] 

[Description of the Related Art] Components to constitute a 

system as a design object are expressed by variable vectors 

consisting of n pieces of elements as follows. 

[0003] 

[ Formula 1 ] 

X v = [X, XJ T (1) 

And the amount to evaluate the characteristics of the 
object is expressed by an evaluation function vector 
consisting of m pieces of elements as follows. 
[0004] 
[Formula 2] 

F v (Xv) = [F.CXJ F n (Xv)] T (2) 

[0005] 
[Formula 3] 

F*(Xv> = w k {f k (X v )-f k , tar } (3) 

Where, suffix v is added to indicate Xv and Fv are 
vectors (or matrix). Further, fk(Xv) , fk and tar are the 
values of the evaluation function to indicate respective 
characteristic values, and a target value thereof, and wk is 
a weighting coefficient. A constraint value may be included 
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in these evaluation functions in addition to the 
characteristic value intended for a desired value. In the 
attenuation least square method, a solution is generally 
obtained so that the amount of deviation (error) from a 
target value of each evaluation function value is minimized, 
and an objective function expressed by the following formula 
is used for a single evaluation scale. 
[0006] 
[Formula 4] 

(MXJ = F V T (X V )F V (X V ) (4) 

The objective function expressed by the formula (4) is 
referred to as a merit function in the field of the lens 
design, in particular. By applying the Newton -Paphson 
method to the formula (4), the necessary information of Xv 
for a local minimum value is given as follows. 
[0007] 
[Formula 5] 

(A v T A v +i; f i H vi ) AX v =-A v t Fv0 (5) 

Where, Av is a Jacobian (Jacob!) matrix of Fv(Xv) , and 
defined by Formula (6). 
[0008] 
[ Formula 6 ] 
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a xi 


d X2 




d F2 


a F2 


Av= 


d xi 


d X2 










d xi 


a X2 



(6) 



a Fi a Fi aFi 

a X n 

a F 2 

a X n 
d Fm 

a xn 

Hvi is the Hesse matrix of the i-th evaluation function, 
and defined by Formula (7). 
[0009] 
[Formula 7] 



Hvi = 



a 2 Fi 


d*Fi 


a*Fi 


a 2 xi 


a xi ax2 


a xi a xn 


a*Fi 




a 3 Fi 


d X2 dXl 


d X2 5X2 


a X2 #Xn 


a 2 Fi 


• 

d z F\ 


■ 


a Xn 5X1 


6 Xn 3X2 





(7) 



In Formulae (6) and (7), Fi means an i-th element 
Fi(Xv) of the evaluation function vector Fv(Xv) (Formula 
(2)). AXv is a solution vector to indicate fluctuation of 
the variable vector to the next step, and expressed as 
follows . 

Generally, in the descent method, an optimum solution 
is obtained by solving Formula (10) by using a diagonal 
matrix Dv (Formula (9)) with a parameter (a non-linear 
correction term) to correct non-linear components of the 
secondary or higher differential as elements thereof in 
place of directly solving Formula (5). 
[0010] 
[Formula 8] 
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AXy = [AX, AXJ T 

[0011] 
[Formula 9] 



(8) 



Dv= 



di 0 



(9) 



[0012] 

[Formula 10] 



(A/Ay + D v )AXv = -Av T F v 



(10) 



In the actual numerical calculation, each element a i:J of 
Av is calculated not by using a strict differential 
coefficient value (Formula 11) shown by Formula (6), but by 
using the difference value (Formula 12). 
[0013] 

[Formula 11] 

a„ = SFi/axj (11) 

a 13 = bF ± /b*i (12) 
[0014] 

[Problems to be Solved by the Invention] However, the value 
of the optimum variable difference amount 6xj is diversified 
according to the characteristic of the object system, 
raising a problem in that, if an inadequate value is given, 
errors occur in the numerical computation, many iterative 
calculations are required, and the optimum solution cannot 
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be efficiently obtained. 

[0015] The present invention is achieved in light of this 
point, and an object of the present invention is to provide 
a design supporting device and a design supporting method 
capable of obtaining the optimum solution in a constantly 
efficient manner by adequately setting the variable 
difference amount. 
[0016] 

[Means for Solving the Problems] In order to achieve the 
above-described object of the present invention, a design 
supporting device to support the design of a lens or the 
like by a descent method for locally changing and optimizing 
a plurality of components so that a system of a design 
object comprising the plurality of components has a desired 
characteristic comprises an eigen value distribution 
information operating means to operate eigen value 
distribution information on the product matrix (Av T Av) 
constituted of the product of the Jacobian matrix (Av) with 
the variable amount of the characteristic of a system to the 
difference amount of the variables corresponding to the 
plurality of components as the elements thereof and the 
transposed matrix (Av T ) of the Jacobian matrix, and a 
difference amount setting means to set the difference amount 
of the variables by using the eigen value distribution 
information. 
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[0017] In order to achieve the same object of the present 
invention, in a design supporting method to support the 
design of a lens or the like by a descent method for locally 
changing and optimizing a plurality of components so that a 
system of a design object comprising the plurality of 
components has a desired characteristic, eigen value 
distribution information is operated on the product matrix 
( Av T Av ) constituted of the product of the Jacobian matrix 
(Av) with the variable amount of the characteristic of the 
system to the difference amount of the variables 
corresponding to the plurality of components as the elements 
thereof and the transposed matrix (Av T ) of the Jacobian 
matrix, and the difference amount of the variables is set by 
using the eigen value distribution information. 
[0018] Further, preferably in the design supporting device 
or the design supporting method, the difference amount of 
the variables is set so that the initial value of the non- 
linear correction term corresponding to the variables to an 
optimizing objective function is present between the minimum 
value and the maximum value of the eigen value distribution 
of the product matrix, the difference amount of the 
variables is set so that the median of the eigen value 
distribution of the product matrix is substantially matched 
with the median of the initial value of the non- linear 
correction term corresponding to the variables to an 
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optimizing objective function, the difference amount of the 
variables is set so that the mean of the eigen value 
distribution of the product matrix is substantially matched 
with the median of the initial value of the non-linear 
correction term corresponding to the variables to an 
optimizing objective function, the difference amount of the 
variables is set so that the mean of the eigen value 
distribution of the product matrix is substantially matched 
with the median of the initial value of the non- linear 
correction term corresponding to the variables to an 
optimizing objective function, the difference amount of the 
variables is set so that the median of the eigen value 
distribution of the product matrix is substantially matched 
with the mean of the initial value of the non- linear 
correction term corresponding to the variables to an 
optimizing objective function, the difference amount of the 
variables is set so that the mean of the eigen value 
distribution of the product matrix is substantially matched 
with the mean of the initial value of the non-linear 
correction term corresponding to the variables to an 
optimizing objective function, and the difference amount of 
the variables is set by using eigen value distribution 
information on the product matrix and predetermined 
weighting information set corresponding to the eigen value 
distribution. 
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[0019] In addition, in the design supporting device or the 
design supporting method, the "square value distribution of 
the singular values of the Jacobian matrix (Av)" may be used 
in place of the "eigen value distribution of the product 
matrix (Av T Av) " . 
[0020] 

[Operation] According to the device of Claim 1 or according 
to the method of Claim 15, eigen value distribution 
information on the product matrix (Av T Av) of the Jacobian 
matrix (Av) with the variable amount of the characteristic 
of the system to the difference amount of the variables 
corresponding to a plurality of components as the elements 
thereof and the transposed matrix (Av T ) of the Jacobian 
matrix, and the difference amount of the variables is set by 
using the eigen value distribution information. 
[0021] According to the device of Claim 8 or according to 
the method of Claim 22, square value distribution 
information on the singular values of the Jacobian matrix 
(Av) with the variable amount of the characteristic of the 
system to the difference amount of the variables 
corresponding to a plurality of components as the elements 
thereof is operated, and the difference amount of the 
variables is set by using the square value distribution 
information. 
[0022] 



- 20 - 



[Embodiments] Embodiments of the present invention will be 
described with reference to the drawings. 

[0023] Fig. 1 shows a configuration of a design supporting 
device according to an embodiment of the present invention. 
This device comprises a CPU (a Central Processing Unit) 1 to 
perform various kinds of operations to optimize a variable 
corresponding to a component of a system for the design 
object, an input device 2 to input data and operation 
commands to the CPU 1 by a user, a storage device 3 to store 
programs executed by the CPU 1 and data or the like under 
operations, and an output device 4 to display the result of 
operations on a display device or to print the result of 
operations by a printer. These components 1 to 4 are 
connected to each other by a bus 5 . 

[0024] Next, the outline of the processing executed by the 
CPU 1 will be described. 

[0025] Firstly, the product matrix (Av T Av) of the above- 
described Jacobian matrix Av and a transposed matrix Av T 
thereof is subjected to the eigen value break- down as 
follows . 
[0026] 

[Formula 12] 

(A/Av) - V V S V V V T (13) 

Where, Vv denotes an orthogonal matrix, and Sv denotes 
an diagonal matrix of n, n-type with n-eigen values of the 
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product matrix (Av T Av) as elements as shown in Formula (14). 
[0027] 

[Formula 13] 

"Si 0" 

Sv= (14) 

0 s* 

The distribution of n- eigen values si, ... sn which 
are thus obtained is generally changed as the difference 
amount Sxj of the above-described variable is changed. In 
this case, the following method is employed in the present 
embodiment to directly control n-eigen value distribution by 
the difference amount 6x j . 

[0028] Firstly, attention is paid to a linear approximation 
in a vicinity of XvO of the above- described evaluation 
function vector Fv(Xv) as expressed below. 
[0029] 

[Formula 14] 

F V (X V ) = F V (X V0 ) + A v AXv (15) 

Then, the i-th element of the evaluation function 
vector Fv(Xv) will be expressed as follows (Formula (16)). 
[0030] 

[Formula 15] 



F 4 <X T ) =Fi<X vQ )+ (" 



Hi 

AX 



)AX A + • 



+ (- 



5 F { 



L )AX a 



AX. 
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In Formula (16), 6jFi denotes the variable amount by the 
i-th variable of the i-th evaluation function. From the 
relationship of Formula (16), it is proved that completely 
similar optimum operation can be performed even by using the 
matrix Av' and the vector AXv 1 such as in Formulae (17) and 
(18) in place of the Jacobian matrix Av (Formula (6)) and 
the solution vector AXv (Formula (8)). 
[0031] 

[Formula 16] 



Av' = 



6 1F1 5 2F1 ••• <5nFi 
6 1F2 6 2F2 • • • 6 nF2 



S lFm S 2Fix 



5 nFn 



(17) 



[(^)---(^r)]' 



(18) 



By this method, the value of each element of the 
Jacobian matrix can be directly changed by the difference 
amount 6x3 of the variable, and n-eigen values of si to sn 
can be freely controlled. 

[0032] As described above, it is determined whether or not 
the set variable difference amount 5xj is an adequate value 
from the relationship between distribution of n-eigen values 
si to sn by the difference amount 6x3 of the variable which 
is once given and the initial value of the damping factor p, 
and the optimum difference amount is set by changing the 
variable difference amount 6xj so as to obtain a more 
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adequate value by using the information, and the optimal 
design value is efficiently obtained. 

[0033] Next, the present invention will be described more 

specifically by an example of the lens design. 

[0034] In the following specific example, the non-linear 

correction term di is expressed by a common d for simplicity, 

and the diagonal matrix Dv with the non-linear correction 

terms as elements thereof is expressed by Formula (19). 

[0035] 

[Formula 17] 

"d 

d 

Dv = . (19) 

Fig. 2 shows a thin-walled lens model as the design 
object according to the first embodiment of the present 
invention. In Fig. 2, the lens has a thickness so as to 
easily show the embodiment. However, the lens thickness and 
the spacing in the calculation are set to be zero. The 
radius of curvature r3 is used so that the focal distance of 
the entire system is correctly 1 to this lens system, and 
remaining two radii of curvature rl and r2 are defined as 
variables. The tertiary spherical difference coefficient 
and the coma aberration coefficient are employed as two 
evaluation functions. The initial shape is expressed by the 
coordinates of (1/rl, l/r2) as follows. 



- 24 - 



1/rl = -1.8, and l/r2 =0.5 (20) 

The difference amounts 5(l/rl) and 6(l/r2) of these two 
variables (1/rl) and (l/r2) are set in a range specified as 
follows . 
[0036] 

[Formula 18] 

SU/rJ = 1 x 10" 5 ~1 x 10°, 5(l/r 2 ) = 1 x 10" 5 ~1 x 10° 

(21) 

In addition, the variable difference amount to be set 
may be selected from four cases for each variable, i.e. , 
sixteen in total at every 10 2 in the range of Formula (21). 
Generally, in the lens design, the range of the variable 
difference amount to be roughly set is empirically 
understood according to the kind of the variables. 
[0036] In this condition, the Jacobian matrix Av for two 
variables and two evaluation functions in each variable 
difference amount set value is obtained. Two eigen values 
si and s2 of the product matrix (Av T Av) are calculated as 
follows (entirely expressed by Formula (22)). 
[0038] 

[Formula 19] 

6(l/rJ = 1 x 10~ 6 , 5(l/r 2 ) = 1 x 10~ 6 

s x = 8.6x10 13 , s 2 = 2.2xl0 -8 (Fig. 3(a)) 
6(l/rJ - lxlO" 6 , 6(l/r 2 ) = lxlO" 4 

s 1 m 8.0xl0" 10 , s 2 = 2.3xl0" 7 (Fig. 3(b)) 



- 25 - 



eu/rj = 1 x lCT 6 , 6(l/r 2 ) = 1 x 10" 2 

Sj = 9.0 x 10' 10 , s 2 = 2.1 x 10" (Fig. 3(c) ) 
SU/rJ = 1 x 10"\ 6(l/r 2 ) = 1 x 10° 

Si = 5.5 x 10" 9 , s 2 = 6.8 x 10° (Fig. 3(d) ) 
6(l/rJ = 1 x 10" 4 , 6(l/r 2 ) = 1 x 10" 6 

Si = 8.6 x 10" 13 , s 2 = 2.2 x 10" 4 (Fig. 4(a)) 
6(l/rJ = 1 x 10' 4 , 6(l/r 2 ) = 1 x 10" 4 

Si = 8.6 x 10*\ s 2 = 2.2 x 10' 4 (Fig. 4(b)) 
6(l/rJ = 1 x 10"\ 6(l/r 2 ) = 1 x 10" 2 

Si = 8.2 x 10~ 6 , s 2 - 2.3 x 10° (Fig. 4(c)) 
6(1/3^) = 1 x 10" 4 , 6(l/r 2 ) = 1 x 10° 

Si = 5.5 x 10" 6 , s 2 = 6.8 x 10° (Fig. 4(d) ) 
5(l/ri) = 1 x 10" 2 , 6(l/r 2 ) = 1 x 10" 6 

Sj = 8.6 x 10" 13 , s 2 = 2.2 x 10° (Fig. 5(a)) 
6(l/rJ = 1 x 10' 2 , 6(l/r 2 ) = 1 x 10" 4 

Sj = 8.6 x 10" 9 , s 2 = 2.2 x 10° (Fig. 5(b)) 
8Cl/r 1 ) = 1 x 10' 2 . 8(l/r 2 ) = 1 x 10' 2 

s 1 = 2.2 x 10°, s 2 = 8.7 x 10" 5 (Fig. 5(c)) 
6(l/r x ) = 1 x 10"\ 5(l/r 2 ) = 1 x 10° 

s x = 4.3 x 10" 1 , s 2 = 8.5 x 10° (Fig. 5(d)) 
Sd/rJ = 1 x 10°, 5(l/r 2 ) = 1 x 10' 6 

s 2 = 5.8 x 10" 13 , s 2 - 5.2 x 10 +3 (Fig. 6(a) ) 
5(l/r 2 ) = 1 x 10°. 5(l/r 2 ) = 1 x 10" 4 

Si = 5.8 x 10" 9 , s 2 = 5.2 x 10 +3 (Fig. 6(b)) 
SU/rJ = 1 x 10°, 6(l/r 2 ) = 1 x 10" 2 
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s 1 = 5.9 x 1CT 6 , s 2 = 5.2 x 10 +3 (Fig. 6(c) ) 
6(l/rJ = 1 x 10°, 5(l/r 2 ) = 1 x 10° 

s x = 1.5 x 10°, s 2 - 5.2 x 10 +3 (Fig. 6(d) ) 

(22) 

As described above, the eigen value distribution of 
each product matrix (Av T Av) is indicated in Fig. 3(a) to Fig. 
6(d). These figures indicate the eigen value distribution 
operated by the CPU 1 and printed by the output device 4 . 
In each figure, the numerals are exponentially indicated on 
the horizontal axis, and the position of each eigen value 
(si, s2) is indicated by *. Further, the position of the 
initial value dO of the non- linear correction term dO is 
indicated by X on the horizontal axis, and the search range 
of the non-linear correction term dO therearound 
(hereinafter, referred to as "d search range") is indicated 
by [## ... ##]. In the present embodiment, dO is set as 
follows . 

dO = 1.0 (23) 

In addition, the p search range is set as follows. 
[0039] 

[Formula 20] 

dO x 10~ 6 - dO x 0.5 x 10 8 (24) 

In the present invention, an adequate value of the 
difference amount of the variable is set from eigen value 
distribution information on the product matrix (Av T Av) , and 
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it is empirically understood in the field of the lens design 
to set one including the eigen values si and s2 in the d 
search range as many as possible. In this case, as clearly 
shown in Figs. 3 to 6, one including no eigen values si and 
s2 in the d search range is a case shown below (Formula 
(25)). 
[0040] 

[Formula 21] 

■ftd/rj = 1 x 10" 6 , 6(l/r 2 ) = 1 x IO" 6 (Fig. 3(a)) 
Sfl/rJ = 1 x IO" 6 , 6(l/r 2 ) = 1 x IO" 4 (Fig. 3(b)) 

(25) 

Further, one including one eigen value (si or s2) is 
the case shown below (Formula (26)). 
[0041] 

[Formula 22] 
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(26) 

Further, one including both the eigen values si and s2 
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is the case shown below (Formula (27)). 
[0042] 

[Formula 23] 
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(27) 

By setting any one of the variable difference amounts 
shown in Formula (27), iterative calculation to search the 
value of the non-linear correction term d can be efficiently 
performed on the average. 

[0043] In order to set an adequate value of the variable 
difference amount more clearly, one in which the initial 
value dO of the non- linear correction term d is smaller than 
the maximum value and larger than the minimum value of the 
eigen value distribution of the product matrix (Av T Av) . 
[0044] 

[Formula 24] 

6(l/rJ = 1 x IO' 6 , 6(l/r 2 ) 
6(l/rJ - 1 x IO' 4 , 5(l/r 2 ) 
6(l/rJ = 1 x IO" 2 , 5(l/r 2 ) 
Sd/rJ = 1 x IO' 2 , 6(l/r 2 ) 



= 1 x 10° (Fig. 3(d) ) 

= 1 x 10° (Fig. 4(d) ) 

= 1 x IO" 6 (Fig. 5(a) ) 

= 1 x IO* 4 (Fig. 5(b) ) 
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6(l/rJ = 1 x 1CT 2 , 6(l/r 2 ) = 1 x 10° (Fig. 6(c)) 

(28) 

In the present embodiment, this requirement is 
satisfied. 

[0045] In addition, an adequate value of the variable 
difference amount is quantitatively set, one may be set, in 
which the median or the mean of the initial value of the 
non- linear correction term corresponding to a plurality of 
variables rl and r2 is substantially matched with the median 
or the mean of the eigen value distribution of the product 
matrix (Av T Av) . Substantial matching is different by the 
object, and it is empirically understood that matching by 10 2 
to 10 3 is sufficient for the case of the lens design. Since 
the initial value dO of the non- linear correction term d is 
commonly given to every above-described variable, the median 
or the mean of the initial value of the non-linear 
correction term can be expressed as do. In other words, 
when the median or the mean of the eigen value distribution 
of the product matrix (Av T Av) is expressed by sm, the 
difference amount 5(l/rl) and 6(l/r2) of the variable set at 
that time if^O satisfies the relationship shown below. 
[0046] 

[Formula 25] 

s m x 10" 3 ~s m x 10" 2 2S^»^s s m x 10 2 ~sm x 10 3 (29) 

In the present embodiment, there are two variables, and 
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the median and the mean of the eigen value distribution of 
the product matrix (Av T Av) are matched with each other. When 
checking the mean of the eigen value distribution 
corresponding to the variable difference amount, it is 
proved that the relationship of Formula (29) is satisfied if 
being set like Formula (28). 

[0047] To set the variable difference amount more 
automatically, weighting as shown in Fig. 7 may be performed 
so that the evaluation value becomes higher as more eigen 
values si and s2 are included in the d search range. In Fig. 
7 , the value on the horizontal axis indicates each eigen 
value, and the value on the vertical axis indicates the 
weight on each eigen value. If the eigen value is included 
in the d search range, 1 is set, otherwise 0 is set. The 
variable difference amount may be set so that the evaluation 
value is much higher thereby. 

[0048] Further, if the variable difference amounts 6(l/rl) 
and 6(l/r2) are set based on the relationship between the 
initial value dO of the non-linear correction term and the 
eigen value distribution, continuous weighting may be 
performed so that the value becomes maximum at the initial 
value dO of the non-linear correction term shown in Fig. 8, 
and the value is gradually decreased as being farther 
therefrom. In addition, weighting information (the 
weighting function) shown in Figs. 7 and 8 is stored in the 
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storage device 3. 

[0049] The result of the optimum operation by setting the 
above-described variable difference amount will be described 
below. 

[0050] Firstly, the number of the eigen values in the d 

search range, the mean of the merit function obtained by the 

result of calculation, and the mean of the iteration number 

required in the calculation will be indicated below. 

[0051] The merit function value is normalized by the value 

at the starting point. 

[0052] 

Number of eigen Mean of merit Mean of iteration 
value function number 

0 0.713 13 

1 0.036 10 

2 0.034 6.7 



As clearly shown from this result, the more the number 
of the eigen values included in the d search range is, the 
smaller the obtained merit function value is, and the 
smaller the iteration number is required. Therefore, the 
design values of the optimum rl and r2 can be obtained by 
setting the variable difference amount so that two eigen 
values are included in the d search range. 

[0053] Next, the result will be shown below for the case in 
which the Initial value dO of the non-linear correction term 
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is present between the minimum value and the maximum value 
of the eigen value distribution and the case in which the 
initial value dO of the non-linear correction term is not 
present therebetween . 
[0054] 

Mean of merit Mean of iteration 
function number 

Present: 0.073 6.6 

Not present: 0.143 10.3 

As clearly shown from this result, by setting the 
variable difference amount so that the initial value dO of 
the non-linear correction term is present between the 
minimum value and the maximum value of the eigen value 
distribution, the obtained merit function value becomes 
small, and the iteration number also becomes small, and the 
optimal design can be efficiently performed. 

[0055] In addition, the result will be shown below for the 
case in which the difference between the initial value dO of 
the non-linear correction term and the mean of the eigen 
value distribution is expressed by the exponent difference, 
and the difference is 10°, 10 3 , 10 4 , 10 7 and 10 8 , respectively. 
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[0054] 



Exponent 
difference 



Mean of merit 
function 



Mean of iteration 
number 



0 



0.066 



6.6 



3 



0.005 



9.0 



4 



0.033 



13.0 



7 



0.469 



13.0 



8 



0.958 



13. 0 



As clearly shown from this result, the merit function 
value is minimized when the exponent difference is about 3, 
and the iteration number becomes less as the exponent 
difference is smaller. Therefore, the optimal design can be 
efficiently performed by setting the variable difference 
amount so that the exponent difference is within a range of 
2 and 3 . 

[0057] Fig. 9 shows a configuration of a lens system as a 
design object according to the second embodiment of the 
present invention . 

[0058] In the present embodiment, twenty two variables 
consisting of lens components are used, and the beam 
aberration is mainly used for the evaluation function to 
indicate the characteristic of the lens. The break-down of 
these twenty- two variables includes eleven surface radii r 
of curvature, and eleven kinds of spacing g, and the 
variable difference amounts are set to an optimum value in 
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terms of the inverse 1/r of the radius of curvature in a 

range as follows. 

[0059] 

[Formula 26] 

5(l/r) = 1 x 10~ 6 ~1 x 10° (30) 

The variable difference amounts are set to an optimum 
value in terms of the spacing g in a range as follows. 
[0060] 

[Formula 27] 

6(g) = 1 x 10- 3 ~1 x 10° (31) 

When each variable difference amount is considered at 
the interval of 10 1 , all possible combinations will be given 
as follows. 
[0061] 

[Formula 28] 

7 11 x 4 11 m 8.3 x 10 15 cases 

Even if all these combinations are checked in 10' 6 
second per combination, it would take an immense amount of 
time of about 263 years. This is one of the most difficult 
problems (a so-called NP- completion type problem) to obtain 
a strict solution in the combination-optimization problem. 
A genetic algorithm (Genetic Algorithm) has been extensively 
known as a method to efficiently solve the approximate 
solution of this kind of problem. The genetic algorithm was 
devised by J. H. Holland in 1960s, and this is a method for 
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efficiently solving an optimization problem by simulating 
the essence of the biological evolution as it is on the 
computer. 

[0062] More specifically, this method is executed by the 
procedure as shown in Fig. 12. Firstly, in Step SI, the 
problem is modeled, and the gene type is set. For example, 
a numerical string (10100110) is used for the gene type. 
Next, a large number of individuals of different genes 
(numerical strings in this case) are formed (for example, 
10110001, 01101010, 11101100 ...) (Step S2). And, the 
adaptability of the formed individuals is evaluated (Step 
S3). In other words, each individual (series) is evaluated 
by using the predetermined evaluation function, and the 
individuals (series) of low evaluation will be selected in 
Step S4. 

[0063] In Step S5, the not-selected individuals are 
multiplied to compensate the number of selected individuals, 
and the crossing-over (Step S6) and the mutation (Step S7) 
are performed. If the adaptability of the obtained 
individuals is insufficient. Steps S3 to S7 will be repeated. 
Here, the crossing- over is performed by switching specified 
parts by selecting specified gene pairs, and the mutation is 
performed by changing some parts of the gene at some 
probability. 

[0064] When the above -described processing is applied to 
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the setting of the optimum variable difference amount, a 
point is how the gene type is set and how the adaptability 
is evaluated. 

[0065] In the present embodiment, firstly, the series in 
which an arbitrary pair of difference amount set values are 
expressed by the exponent are employed. In other words, the 
first to eleventh number of the series from the left to the 
right are set as the exponent value of the difference amount 
on each radius r of curvature, the twelfth to twenty- second 
number of the series are set to be the exponent value of the 
difference amount on each spacing g, all the difference 
amounts of 1/r are set to be 1 x 10" 6 , and all the difference 
amounts of g are set to be 1 x 10" 3 , the gene row to indicate 
this state is expressed as follows. 

-6-6-6-6-6-6-6-6-6-6-6-3-3-3-3-3-3-3-3-3-3-3 (32) 
[0066] On the other hand, information on the eigen value 
distribution of the product matrix (Av T Av) is used for the 
evaluation of the adaptability. More specifically, the 
result of higher evaluation will be obtained when the 
initial value dO of the non-linear correction term dO is 
smaller than the maximum of the eigen values of the product 
matrix (Av T Av) , and larger than the minimum thereof. 
[0067] When the adaptability is evaluated more strictly, 
the median or the mean of the initial value dO is compared 
with the median or the mean of the eigen value distribution. 
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and higher evaluation will be obtained when the initial 
value is substantially matched with the median or the mean 
of the eigen value distribution. 

[0068] In the present embodiment, weighting of the 
distribution as shown in Fig. 10 is performed, and the 
evaluation value becomes higher when the eigen value is not 
outside the d search range but within the range, and higher 
evaluation is obtained when the eigen value is close to the 
initial value dO of the non- linear correction term though it 
is within the d search range. 

[0069] Description will be given by a specifically 
numerical example. Firstly, in the variable difference 
amount , the radius r of curvature and the spacing g of the 
plane are set to be common as follows. 
[0070] 

[Formula 29 ] 

6(l/r) m 1 x 10" 6 

5(g) = 1 x 10" 3 (33) 

In this condition, the eigen value distribution is 
examined as shown in Fig. 11(a). This state forms a 
starting point, and the weighting distribution in Fig. 10 is 
multiplied by the eigen value distribution in Fig. 11(a) to 
obtain the evaluation value, and, processing is performed to 
set an adequate variable difference amount by the genetic 
algorithm. The variable difference amount set value 
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obtained as the result is as shown below for the radius r of 

curvature and the spacing g, respectively. 

[0071] 

[Formula 30] 

6(l/r) - 1 x lO" 5 -! x 10* 2 

5(g) = 1 x 10* 3 ~1 x 10- 1 (34) 

The eigen value distribution corresponding to the 

variable difference amount set value (34) is shown in Fig. 

11(b). As clearly shown from the comparison between Fig. 

11(a) and Fig. 11(b), it is shown that the eigen value 

distribution is gathered around the initial value dO of the 

non-linear correction term dO (= 1.0) by the above-described 

processing . 

[0072] Comparison of the normalized merit function value 
and the iteration number required for the calculation before 
and after the above -described processing is as follows. 
[0073] 

Merit function Iteration number of 
value calculation 

Before processing 5.0 x 10" 3 3 

After processing 1.3 x 10° 3 

As clearly shown from this result , according to the 

present embodiment, the convergence efficiency of the 

descent method can be enhanced. 

[0074] Next, the outline of the processing executed by the 
CPU 1 in the third embodiment of the present invention will 
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be described. 

[0075] Firstly, the Jacobian matrix Av is subjected to the 

eigen value break-down. 

[0076] 

[Formula 31] 

^ = U V P V V V T (35) 

Where, UvVv denote orthonormal matrixes of m, m type 
and n, n type, respectively, and when m > n, Pv is expressed 
as follows. 
[0077] 

[Formula 32] 
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pi 
o 



Pn 



(36) 



When m ^ n, Pv is expressed as follows, 
[0078] 

[Formula 33] 

pi 0 



0* 



Pm 



(37) 



Pv is a singular value matrix with n- or m- singular 
values of the Jacobian matrix Av as the elements thereof. 
The product Pv T Pv of the singular value matrix is expressed 
as follows. 
[0079] 
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[Formula 34] 

(38) 

The distribution of the square values pi 2 , ... pn 2 of 
the thus -obtained n- singular values is changed by the 
difference amount 6xj of the above -described variables. In 
this condition, in order to directly control the 
distribution of the square values pi 2 , . . . pn 2 of the thus- 
obtained n- singular values by the difference amount 5x j , a 
method (Formulae (15) to (18)) similar to that of the first 
embodiment is employed also in the present embodiment. 
Therefore, the square values pi 2 , ... pn 2 of the thus- 
obtained n- singular values can be freely controlled by the 
difference amount 5x j . 

[0080] As described above, it is determined whether or not 
the set variable difference amount is adequate from the 
relationship of the square value distribution of the 
distribution of n- singular values of the once-given 
variables by the difference amount 5xj and the initial value 
of the non- linear correction term d, and optimum difference 
amount is set by successively changing the variable 
difference amount Sxj so as to obtain a further adequate 
value by using the information so that the optimal design 
value is efficiently obtained. 



p v I Pv= 



Pi 



Pn* 
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[0081] The above -described method is applied to the same 
lens system identical to that of the first embodiment, and 
the square values pi 2 and p2 2 of the singular values of the 
Jacobian matrix Av corresponding to Formula (22) are 
calculated as follows (the entire values are expressed by 
Formula ( 39 ) ) . 
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As clearly shown above, the distribution of the square 
values of the singular values of the Jacobian matrix Av is 
totally same as the distribution of the eigen values of the 
product matrix ( Av T Av ) . In other words, the square values 
pi 2 and p2 2 of the singular values of the Jacobian matrix Av 
are equal to the eigen values of the product matrix (Av T Av) . 
[0083] Therefore, in the first embodiment, replacement of 
the singular values of the eigen values si and s2 of the 
product matrix (Av T Av) by the square values pi 2 and p2 2 
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corresponds to the present embodiment, and similar effect to 
the first embodiment can also be demonstrated by the present 
embodiment . 

[0084] Further, similar replacement can also be performed 

even in the second embodiment . 

[0085] 

[Advantages] As described above in detail, according to the 
design supporting device of Claim 1 , or according to the 
design supporting method of Claim 15, eigen value 
distribution information on the product matrix (Av T Av) 
constituted of the product of the Jacobian matrix (Av) with 
the variable amount of the characteristic of the system to 
the difference amount of the variables corresponding to a 
plurality of components as the elements thereof and the 
transposed matrix (Av T ) of the Jacobian matrix is operated, 
and the difference amount of the variables is set by using 
the eigen value distribution information. Therefore, the 
optimum value by the descent method is efficiently operated, 
and the optimal design value can be efficiently obtained. 
[0086] Further, according to the design supporting device 
of Claim 8 or according to the design supporting method of 
Claim 22, square value distribution information on the 
singular values of the Jacobian matrix (Av) with the 
variable amount of the characteristic of the system to the 
difference amount of the variables corresponding to a 
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plurality of components as the elements thereof, and the 
difference amount of the variables is set by using the 
square value distribution information, and similar effect 
can be demonstrated thereby. 
[Brief Description of the Drawings] 

[Fig. 1] Fig. 1 shows a configuration of a design 
supporting device according to an embodiment of the present 
invention. 

[Fig. 2] Fig. 2 shows a configuration of a first lens 
system as a design object. 

[Fig. 3] Fig. 3 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) of a 
product matrix of the Jacobian matrix and a transposed 
matrix thereof . 

[Fig. 4] Fig. 4 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) of a 
product matrix of the Jacobian matrix and a transposed 
matrix thereof . 

[Fig. 5] Fig. 5 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) of a 
product matrix of the Jacobian matrix and a transposed 
matrix thereof . 

[Fig. 6] Fig. 6 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) a 
product matrix of the Jacobian matrix and a transposed 
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matrix thereof . 

[Fig. 7] Fig. 7 shows an example of weighting information 
(a weighting function). 

[Fig. 8] Fig. 8 shows an example of weighting information 
(a weighting function) . 

[Fig. 9] Fig. 9 shows a configuration of a second lens 
system as a design object. 

[Fig. 10] Fig. 10 shows an example of weighting information 
(a weighting function). 

[Fig. 11] Fig. 11 shows distribution of the eigen value (a 
square value of a singular value of a Jacobian matrix) a 
product matrix of the Jacobian matrix and a transposed 
matrix thereof. 

[Fig. 12] Fig. 12 shows a flowchart to indicate the 
procedure of a genetic algorithm. 

[Reference Numerals] 

1 CPU 

2 input device 

3 storage device 

4 output device 
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FIG. 1 

2 INPUT DEVICE 

3 STORAGE DEVICE 

4 OUTPUT DEVICE 

FIG. 7 

1 WEIGHTING VALUE 

2 WEIGHTING DISTRIBUTION 

3 NUMERICAL VALUE 

4 dO VALUE 

5 d SEARCH RANGE 

FIG. 8 

1 WEIGHTING VALUE 

2 WEIGHTING DISTRIBUTION 

3 NUMERICAL VALUE 

4 d SEARCH RANGE 

FIG. 10 

1 WEIGHTING VALUE 

2 WEIGHTING DISTRIBUTION 
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3 NUMERICAL VALUE 

4 dO VALUE 

5 d SEARCH RANGE 

FIG. 12 

51 SET GENE TYPE 

52 SET GROUP 

53 EVALUATE ADAPTABILITY 

54 SELECTION 

55 MULTIPLICATION 

56 CROSSING-OVER 

57 MUTATION 

58 IS ADAPTABILITY SUFFICIENT? 
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